Simultaneous Onset of Condensation of Molecules and Atoms 
in an Attractive Fermi Gas of Atoms 
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The self-consistent equations for the order parameters of Bose-Einstein condensation (BEC) of 
molecules and Bardeen-Cooper-Schrieffer (BCS) condensation of atoms in a Fermi gas of atoms 
with an attractive two-body interaction between atoms have been derived within the Hartree-Fock- 
Bogoliubov approximation from the path integral representation of the grand partition function. 
We have found that the order parameters for BEC and BCS are proportional to each other, which 
implies that BEC and BCS onsets simultaneously. We have also found that the common critical 
temperature of BEC and BCS increases as the average number of molecules increases and that the 
atom-molecule coupling enhances the common critical temperature. 
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The crossover regime [1HSSHE0@ ° f B ° Se " 
Einstein condensation (BEC) of molecules and Bardeen- 
Cooper-Schrieffer (BCS) condensation of atoms in a 
Fermi gas of atoms has recently been an active research 
field since it contains abundant information on the inter- 
play of two kinds of quantum fluids (the Bose and Fermi 
quantum liquids) that exist in the nature. Progress in 
understanding of this crossover regime will certainly fur- 
nish significant implications and shed much light on our 
understanding of a plethora of macroscopic quantum phe- 
nomena including high temperature superconductivity. 
However, despite that great efforts have been put into it, 
the nature of condensation on either side of the crossover 
regime still remains an open question. 

In this Letter, we report our results obtained in our 
study of Bose-Einstein condensation of molecules and 
BCS condensation of atoms on the BCS side in a Fermi 
gas of atoms with an attractive two-body interaction be- 
tween atoms (an attractive Fermi gas of atoms). We 
will first derive a set of self-consistent equations for the 
order parameters of BEC and BCS in the Hartree-Fock- 
Bogoliubov approximation from the path integral repre- 
sentation of the grand partition function. Then, from the 
self-consistent equations we will draw the central conclu- 
sion in this Letter that BEC of molecules and BCS of 
atoms onset simultaneously in an attractive Fermi gas 
of atoms. We also calculate the common critical tem- 
perature of BEC and BCS and discuss its dependence 
on the fraction of atoms bound into molecules, on the 
two-body interaction between atoms, and on the atom- 
molecule coupling. 

In the literature, the prime examples of Fermi gases 
of atoms are the gases of 40 K and 6 Li atoms, with the 
gas of 40 K atoms most intensively studied so far. The 
magnetic-field Feshbach resonance 0, is used to vary 
the sign and strength of the two-body interaction be- 
tween atoms so that the BEC-BCS regime is conveniently 
accessed. In our present work, we are concerned with 
the BCS side in which the two-body interaction between 



atoms is attractive. On the BCS side, there exist an 
open channel and a closed channel, with bound states 
only in the closed channel. The bound-states in this case 
are referred to as unstable (unphysical) bound states. A 
resonance molecule refers to the entity formed by two 
atoms in such a bound state in real space. If the size of a 
bound state is small in real space, the bound state is then 
a conventional molecule. The size of a bound state can 
be varied by varying the strength of the two-body inter- 
action between atoms through the Feshbach resonance. 
As its size increases, the bound state becomes better lo- 
calized in momentum space and it can then be referred 
to as a preformed Cooper pair. 

A Fermi gas of atoms with an attractive two-body in- 
teraction is described by the Hamiltonian 

H = 2Ai/^aj>p + ^(e k -/i)c£ .c fcff 
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where 2Af = 1(y — /i) is the offset energy of a molecule 
relative to the chemical potential /i, a p and a p are cre- 
ation and annihilation operators of a molecule of mo- 
mentum p and energy 2Av, c£ and Cka are creation and 
annihilation opperators of a fermionic atom of momen- 
tum k, spin a, and energy efc = h 2 k 2 /2m, U (< 0) is the 
two-body interaction between atoms, g is the effective 
atom-molecule coupling constant, and N is the average 
number of single atoms. 

Our starting point is the path integral representa- 
tion 0,0] of the grand partition function Z = Tre _/3ff 
with H given in Eq. (JJJ . To obtain the path integral rep- 
resentation of the grand partition function, we introduce 
the complex variables 92* (r) and f p (r) for the creation 
and annihilation operators al(r) and a p (r) of molecules, 
and the Grassman variables 4'ka( T ) an d V'fco-fV) for the 
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creation and annihilation operators cj^r) and Cfc CT (r) 
of atoms. Replacing a p (r), a p (r) by <£> p (t), f P (r) and 

c L( r )' c feo-(T) by V>Lt( t )' V'fecr(' r ) m tne Hamiltonian in 
Eq. so that the Hamiltonian becomes a function of the 
above-introduced complex and Grassman variables and is 
denoted by H(ip,ip), and then substituting H(ip,?p) into 

Z = J D[(p*,ip;tp*,ip] T expj-y dr <P* p (t)-^(p p (t) 

+ V>LM|:Vw(r)-/%,^)]} ) (2) 

we obtain the integral representation of the grand parti- 
tion function. The integration measure in Eq. (J2J is given 
by 
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with Ar = (3/M. 

The complex variable f P (r) introduced above acts as 
the order parameter of Bose-Einstein condensation of 
molecules. To proceed further, we use Fourier trans- 
formations of f p {r) and V-'fco-M given by f p {r) = 
E iUm e~ iUmT <P P and ^(t) = E u /" u " T t [where 
iuj m = ilmm j (5 and = i(2n + l)ir/P, with m and 
n being integers, are Matsubara imaginary frequencies 
for bosons and fermions, respectively, and p = (p,iu> m ), 
k = (k,iu) n )] to rewrite the path integral representation 
in imaginary frequency space. For the purpose of decou- 
pling the two-body interaction between atoms through 
the Hubbard-Stratonovich transformation 0, 0], we 
rewrite it as 
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The field introduced in implementing the Hubbard- 
Stratonovich transformation is denoted by A p and it cor- 
responds to the order parameter of BCS condensation of 
atoms ^2 ■ For discussion of Bose-Einstein condensation 
of molecules and BCS condensation of atoms, we need 
retain only zero bosonic momentum-frequency terms in 
the atom- molecule coupling and in the decoupled terms. 
Then all nonzero momentum-frequency order parameters 
can be integrated out, with only the order parameters of 
zero momentum and frequency [they are now simply de- 
noted by tp and A, respectively] remaining in the path 
integral representation of the grand partition function. 

The next algebraic manipulation to perform is the in- 
tegration over the Grassman variables. To integrate out 



the Grassman variables, we first bring the terms contain- 
ing the Grassman variables as a whole into a diagonal 
form. This manipulation is similar to the diagonaliza- 
tion of the Hamiltonian of a system of fermions . We 
found that it is convenient to introduce a column matrix 
of two Grassman variables and to write the concerned 
terms in a matrix form. Then we diagonalize the matrix 
appearing in the matrix form of the concerned terms. 
The transformation matrix that transforms the to-be- 
diagonalized matrix into a diagonal matrix is composed 
of the orthonormal eigenvectors of the to-be-diagonalized 
matrix. The new Grassman variables in terms of which 
the concerned terms are expressed in a diagonal form can 
be easily obtained by using the transformation matrix. It 
can be directly verified that the Jacobian of the change 
of variables from the old Grassman variables to the new 
Grassman variables is unity. 

The final simplification we can do to the path integral 
representation of the grand partition function is to per- 
form the summations over the imaginary frequencies iuj m 
and iw n . These summations can be easily performed by 
considering proper contour integrals [HEi. With the 
summations over ioj m and ioj n done, we obtain 



Z = N 2 f3\U\ 
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Notice that the prefactor N 2 j3\U\ in Eq. © arises from 
the change of variables <p — > \fNip and A — > \/~N (3\U\A 
before performing the summations over iuj m and iu) n . 
The purpose of this change of variables is to make the 
order parameters ip and A intensive. 

An approximate explicit expression of Z can be ob- 
tained by making use of the steepest-descent method 0] 
in which S is expanded as a Taylor series in ip and A 
about the stationary point of S. As an approximation, 
this series is terminated at the second order. The in- 
tegration over if and A can then be performed and an 
approximate explicit expression of Z is obtained. The 
stationary value of S is obtained by evaluating its value 
at the stationary point determined by the conditions that 
dS/dip* — and dS/dA* = and their complex conju- 
gates. From the stationary conditions, we obtain the 
self-consistent equations for ip and A 
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The self-consistent equations for the complex conjugates 
of ip and A can be obtained simply by taking the complex 
conjugation of the above two equations. 

From Eqs. JSJ, we immediately see that ip and A are 
proportional to each other 



a 



2Av 



A. 



(9) 



The proportionality between tp and A implies that they 
become nonzero at the same temperature as the temper- 
ature is lowered from a value higher than the critical tem- 
perature. In other words, Bose-Einstein condensation of 
molecules and BCS condensation of atoms occur simulta- 
neously. This result is of great significance in that it re- 
veals to us that both molecules and atoms condense in an 
attractive Fermi gas of atoms and furthermore that their 
condensation onsets at a common temperature. Physi- 
cally, the simultaneity of condensation of molecules and 
atoms is due to the presence of the atom-molecule cou- 
pling. 

We now turn to the evaluation of the common criti- 
cal temperature of condensation of molecules and atoms. 
The common critical temperature T c can be determined 
from any one of the two self-consistent equations in 
Eqs. (JHJ in conjunction with the relation in Eq. ©. Ex- 
plicitly, we have 



1 = 
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Notice that, since the chemical potential appears in 
Av = v — /i, the chemical potential \x has to be de- 
termined self-consistently. To determine fi, we set up 
a self-consistent equation for it by using the relation 
2M + N = —dF/dfi in which F is the Helmholz free 
energy and M is the average number of molecules. This 
relation can be easily verified by using F = — /3~ l \n.Z 
and Z = Tre _/3H with H given in Eq. JTJ. With Z ex- 
pressed as Z = e , the above relation can be written as 
2M+N = -/3- 1 dS/d(i. In the Hartree-Fock-Bogoliubov 
approximation, it is satisfactory to approximate S with 
its stationary value. With the utilization of such an ap- 
proximation to S, we obtain the following self-consistent 
equation for the chemical potential /i c at T = T c 



2rja_ \ - 1 
N ^ e^" + 1 ' 

k 



,2/3 c Av _ 1 2(3 c Au 



(11) 

where (3 C = l/ksT c , r\ m is the ratio of the number of 
atoms bound into molecules to the total number of atoms, 
rj m = 2M/(2M + N), and r\ a is the ratio of the number of 
single atoms to the total number of atoms, r\ a = N/ (2M+ 
N). Notice that, in deriving the self-consistent equation 
for jU c , the zero-point contributions must be omitted |20j |. 

Before we proceed to solve for the common critical tem- 
perature T c from Eqs. (|10fl and i|ll|) . we first specify val- 



ues for several parameters. Firstly, an energy cutoff is re- 
quired in Eq. (|10|l . This cutoff indicates the energy range 
of atoms that participate in pairing. We take this energy 
cutoff to be the Fermi energy Ep since all atoms are al- 
lowed to participate in pairing . For a typical density of 



icipatc j 

n = 10^ u m" 3 0,E3, E F « 1.25 (jK for 4U K atoms. Sec- 
ondly, the range of the atom-molecule coupling g we will 
consider is between 0.1-Ef and Ep. Thirdly, for the two- 
body interaction between atoms, we use U = —0.2Ef- 
Lastly, for the molecular offset energy, we use v = Ef- 
The above-quoted values of parameters are in consistency 
with the experiment Q] and with the previous theoretical 
analysis |ll|. 

Shown in Fig. QJa) [the solid line] is the dependence of 
the common critical temperature T c on the fraction r\ m 
of atoms bound into molecules, with the dependence on 
T] m of the chemical potential fi c at T c shown in Fig. ^b) . 
It is seen that, as r\ m increases, T c increases almost lin- 
early, with T c w 0.12E F at rj m = and T c sa 0.63^ at 
Vm = 0.25. We thus conclude that the presence of more 
molecules effectively raises the common critical temper- 
ature of condensation. The dependence of T c on rj m can 
be well fitted with the power law ksTc/Ep = a + brf^ 
[shown in Fig. QJa) as a dashed-line] with a w 0.120, 
b w 1.569, and 7 w 0.815. From Fig. ^b), we can see 
that, as rim increases, the value of the chemical potential 
fi c at T c also increases, with the increase slowing down 
at large values of r\ m . The dashed-line in Fig.QJb) is the 
fit to the power law ksTc/Ep = c + drj m with c = 0.984, 
d = 1.824 x 10~ 2 , and S = 0.233. 
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FIG. 1: The dependence of the common critical temper- 
ature Tc [(a)] and the chemical potential /j, c at T c [(b)] on 
the fraction of atoms bound into molecules for g — O.IEf, 
U = — 0.2Ef, and v — Ef- The dashed lines are power-law 
fits discussed in the text. 

Both the two-body interaction between atoms and the 
atom-molecule coupling have significant influences on the 
common critical temperature T c as shown in Fig. [21 for 
Vm = 0.15 and v = Ep. At these values of rj m and v, the 
values of T c /Ep fall in the vicinity of 0.5Ep for values of 
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\U\/Ep and g/Ep within the range of 0.1 ~ 1.0, which 
is in consistency with the experiment Q. As \U\ or g 
increases, T c increases monotonically. It has been found 
that both the T c -versus-|[/| and T c - versus- g curves can 
be well fitted with a power law fcsT" c / Ep — a + (3x p . 
For the Inversus- curve, x = \U\/E F , a = 0.452, 
13 = 7.85 x 10" 3 , and p = 2.410; for the T c -versus-g 
curve, x = g/E F , a = 0.435, (3 = 0.257, and p = 1.327. 
These fits are shown in Figs. Ela) and Gib) as dashed- 
lines. From Fig. it can be seen that the influence on 
T c of the atom-molecule coupling is stronger than that of 
the two-body interaction between atoms. The reason for 
this behavior is the fact that g appears in the quadratic 
form in the effective coupling constant (g 2 /2Aiy + \U\) 
[the first factor on the right hand side of Eq. QlOfl] , while 
\U\ appears only in the linear form. 
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FIG. 2: The effects of the two-body interaction [(a)] and the 
atom-molecule coupling [(b)] on the common critical temper- 
ature T c . The dashed lines are the power-law fits discussed in 
the text. 

In encapsulation, in terms of the self-consistent equa- 
tions for the order parameters derived from the path 
integral representation of the grand partition function, 
we have studied Bose-Einstein condensation of molecules 
and BCS condensation of atoms in an attractive Fermi 
gas of atoms as realized in experiments 0, 0, El by mak- 
ing use of the magnetic-field resonance. We have found 
that molecules and atoms start to condense at the same 
critical temperature, that the common critical tempera- 
ture increases as the number of molecules increases, and 
that the atom-molecule coupling raises the common crit- 
ical temperature. 

Our findings have a number of important implications 
both to the BEC-BCS crossover in particular and to su- 
perfluidity (superconductivity) in general. The fact that 
both molecules and atoms condense and that they have 
the very same critical temperature resolves the contro- 
versy over the nature of condensation on the BCS side 
of the BEC-BCS crossover regime and manifests that the 
condensate in an attractive Fermi gas of atoms is actu- 



ally composed of a Bose-Einstein condensate and a BCS 
condensate. The increase of the common critical temper- 
ature with the fraction of atoms bound into molecules 
indicates that high temperature (possibly room temper- 
ature) superfluidity (superconductivity) is likely to be 
achieved in a system in which Bose-Einstein condensation 
of bosons and BCS condensation of fcrmions can occur 
simultaneously. The enhancement of the common criti- 
cal temperature by the atom-molecule coupling suggests 
that high temperature superfluidity (superconductivity) 
is likely to be discovered in a system with a strong boson- 
fermion interaction in addition to that both bosons and 
fermions in the system can undergo, respectively, Bose- 
Einstein and BCS condensations. 
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